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Abstract. Beyond the main sequence solar type stars undergo extensive mass loss, providing an 
environment where planet and brown dwarf companions interact with the surrounding material. To 
examine the interaction of substellar mass objects embedded in the stellar wind of an asymptotic 
giant branch (AGB) star, three dimensional hydrodynamical simulations at high resolution have 
been calculated utilizing the FLASH adaptive mesh refinement code. Attention is focused on the 
perturbation of the substellar mass objects on the morphology of the outflowing circumstellar matter. 
In particular, we determine the properties of the resulting spiral density wake as a function of the 
mass, orbital distance, and velocity of the object as well as the wind velocity and its sound velocity. 
Our results suggest that future observations of the spiral pattern may place important constraints on 
the properties of the unseen low mass companion in the outflowing stellar wind. 
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INTRODUCTION 

Evolved giant stars lose a significant amount of mass in the form of stellar winds, 
providing material in regions for which orbiting planets or brown dwarfs interact with 
their immediate surroundings. The mass loss rate of these giant stars, especially during 
the asymptotic giant branch (AGB) phase of the stellar evolution, ranges from 10^^ to 
10^-^ Mq yr^^ with a wind speed of 10-30 km s^^ (Fong et al. 2006). Such rates lead 
to an envelope density at ~ 100 AU distance from the stellar center of lO^^^-lO^^^ 
g cm~^, approaching the typical density of protoplanetary disks (Hayashi 1981; Aikawa 
& Herbst 1999) where planets form and grow by gathering the surrounding matter (e.g., 
Alibert et al. 2005; Hubickyj et al. 2005; Ikoma et al. 2000). The gravitational interaction 
between a protoplanet and its environmental disk has been investigated in great detail. 
Here, a spiral wave forms in the differentially rotating system (Goldreich & Tremaine 
1979; Masset 2008, and references therein), but the corresponding study in the context 
of systems interacting with an outflowing wind around the evolved giant stars is still 
lacking. Although the velocity direction of the background medium is the only difference 
between the rotating and outflowing systems, the entire background structure and the 
interaction aspects markedly differ. In the former case, where the surrounding matter is 
in the form of a differentially rotating disk, the background co-rotates locally with the 
substellar mass component of interest so that the torque at the resonance plays the crucial 



role in producing the spiral. However, in the case of a nearly spherically symmetric 
stellar wind, the motion of the substellar object is perpendicular to the background 
velocity field. The relative velocity of the background gas with respect to the substellar 
object causes a Bondi-Hoyle-Lyttleton (BHL) accretion column to lie behind the orbiting 
body (Bondi & Hoyle 1944; Edgar 2004, and references therein). 

Although starting from a different viewpoint, the accretion flow in a weak gravita- 
tional potential (so-called a density wake) is well formulated by a linear perturbation 
analysis in an initially uniform and static background (Ostriker 1999; Kim & Kim 2007; 
Kim et al. 2008). In particular Kim & Kim (2007) considered an object moving on a 
curvilinear orbit, which leads to a density wake in the form of an Archimedes spiral 
with the opening angle 6 = sin^(cs/Vp), where Vp and Cs represent the orbital and 
sonic speeds. In this contribution we show that the wind from the evolved central star 
pushes the wake outward in the radial direction so as to widen the opening of the spiral 
pattern. This provides the possibility for detecting the spiral pattern at large distances 
from the luminous central star and can be used to probe the presence of the unseen ob- 
jects orbiting around the evolved star, if the density enhancement is sufficient. To predict 
the orbital properties of these substellar components from future observations of the spi- 
ral density wakes, we provide an empirical formula for the density enhancement in the 
spiral pattern. The reader is referred to Kim & Taam (201 1) for a more comprehensive 
discussion. 

SIMULATION SETUP 

Given a continuous stellar mass loss rate M*, a self-consistent configuration of density 
Pw and velocity satisfies the steady hydrodynamic conditions 

V-(PwVw)=0, (1) 

and 

Vw • VV^ = -—VP^- (2) 

Pw 

for the gravitational potential of the central star (Gail & Sedlmayr 1987, Lamers & 
Cassinelli 1999, see also Parker 1958 for the original description for isothermal wind 
solutions). To simplify the problem, we do not consider radiation pressure on dust cou- 
pled with the nearby gas, which is commonly regarded as the wind driving mechanism 
of the evolved stars (e.g.. Winters et al. 2000). Although the driving mechanism affects 
the background wind condition, the response of gas to the gravitational perturbation due 
to the substellar object orbiting around the mass losing star is likely independent of the 
global wind structure as it is determined by the local wind quantities. This has been 
verified by performing two sets of simulations with transonic and supersonic branches 
of Parker's wind solution, confirming that the wake properties described here are deter- 
mined by the wind density pwir), its expansion speed Vw{r), and the sound speed Cs(r) 
on the spot in addition to the orbital and accretion properties of the perturbing object. 

We have investigated 42 hydrodynamic models for a range of object mass M^, size r^., 
orbital speed Vp, distance from the star Vp, as well as stellar mass loss rate M^, wind speed 




FIGURE 1. Comparison of perturbed density distribution between models with the background wind 
speed Vw of (a) 0, (b) 5, and (c) 10 in the units of sonic speed for a background gas characterized by 
Cs — 3 km s^^ The perturbing object is modeled by an object of mass Mp =0.1 Mrr, moving on a circular 
orbit at a distance = 20 AU, marked by a white circle, in the counterclockwise direction with the 
orbital speed Vp = lOcj, currently located at {x,y) = {rp,0). The white line denotes the opening angle of 
the spiral in the vicinity of the orbiting object in a static medium, 9 = sin^' (cs/Vp). The color bar labels 
the perturbed density in logarithmic scale. 



Vw, and sound speed Cj. Specifically, we have verified the dependence of density wake 
on three different velocities: Cs = 3 and 5 km s^^, Vp is in the range of 0.5-lOcs, and Vw 
up to lOcs satisfying the supersonic and transonic branches of isothermal wind models. 
The simulations are performed in three spatial dimensions, but here we only focus on 
the features in the orbital plane. See Kim & Taam (201 1) for the vertical structures. 



SHAPE OF SPIRAL PATTERN 



The outflowing wind affects the shape of the spiral density wake, which is generated 
by the gravitational interaction of the substellar object with the background. As shown 
in Figure 1, the spiral wake pattern opens further with higher wind speeds. To show the 
difference of the opening angle of the spiral, one line is overlaid at the present position of 
the object, showing the angle, G = sin~^ (cs/Vp) , analytically calculated for the model in 
the absence of a wind (Ostriker 1999; Kim & Kim 2007). This expression for the opening 
angle reveals that the pattern propagates in the radial direction with the sonic speed 
as the perturbing object moves laterally with the speed Vp, i.e., the pattern propagation 
speed is Cg (Kim & Kim 2007). Extending this picture to the outflowing models, one can 
glean from Figure 1 the wind speed is related to the pattern propagation speed, which is 
revealed in the opening angle of the spiral arm pattern. 
Simple geometry yields the following relation 

d{r/rp) _ Varm 
d(p Vp ' 

from which the pattern propagation speed Vam relative to the orbital speed Vp of the 
perturbing object can be derived from the shape of the spiral. Numerical differentiation 
at the shock boundary, exhibiting high density, gives the pattern propagation speed of 
Vw + Cs, which is also predicted by linear analysis in the WKB approximation. Here, 
we point out that it is not well described purely by the speed of outflowing matter, 
which is often used by observers (e.g. Mauron & Huggins 2006; Morris et al. 2006), 
especially when the wind speed is comparable to the thermal sound speed. Also it is 
to be noticed that the pattern is not necessarily an Archimedes spiral, in which the arm 
spacing is constant, but its shape possibly changes with the variation of wind speed and 
temperature as a function of distance. 



DENSITY ENHANCEMENT 

To detect the spiral arm pattern around the evolved giant stars, the density enhancement 
(or density jump) of the structure must exceed some threshold value. In our definition 
a = p /pw — 1 represents this perturbed density normalized by the local background 
density pwir), which has the profile close to ~ in most areas. 

Figure 2 exhibits the density enhancement for four different models with different 
combinations of the orbital speed and the wind speed but with a fixed sound speed, as 
well as the mass, radius, and orbital radius of the object, and the stellar mass and mass 
loss rate. We have also checked the effects of all other parameters enumerated here to 
obtain an empirical formula for the minimum value of the density peaks at the shock 
boundary of the spiral arm: 
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FIGURE 2. Density enhancement a of the spiral density wake formed by the gravitational interaction 
of the stellar wind with a substellar object of mass 0.1 M© orbiting around the central star as a function of 
distance normalized by the orbital radius rp = 20 AU. Comparison is achieved between the models with 
the orbital and wind speeds of (a) Vp = 5cs, = lOcg, (b) Vp = lOcs, Vw = lOcs, (c) Vp = 5cs, Vw = 5cs, 
and (d) Vp = lOcs, Vw = 5cs, where the sound speed Cs is 3 km s~' . The lines denote the perturbed density 
profiles along the distance from the star in the direction of the object (solid) and in the opposite direction 
{dashed). See text for the complementary bold solid hne that outlines the peaks. 



where rg = GMp/c^ is the Bondi accretion radius. We mention that this empirical 
formula can be improved by exploring additional parameter space and/or by direct 
analytic approach, although it serves quite well for a large range of parameter space 
as the lower limit of density peaks that could be detected. Note that this equation is well 
checked only for the cases of a supersonic speed of wind, even though the density wakes 
in the transonic wind background appear reasonable in the inner subsonic regime. 

For a fast wind, which is often observed in AGB envelopes, this equation (4) is 
reduced simply to 

r 1 /2 

ocpeak ^ (y ) for Vw > Vp, Cs, (5) 

when the observed arm is located at great distances from the central star relative to the 
orbital radius of the perturbing object. In comparison to the slow wind cases, the density 
peak profiles decrease with distance from the star with a power of —0.5 rather than — 1. 



This is favorable for observers seeking a signature of unseen low mass objects revolving 
about evolved mass losing giant stars. 



DISCUSSION 

Prediction of the Orbital Properties of Planets and Brown Dwarfs 
Deduced from Spiral Structures of AGB Stellar Wind Envelopes 

Equation (5) indicates that for the typical envelope temperature of the cool giants, 
corresponding to the sound speed of 1-4 km s \ the peak density jump for objects of 
1 Mj, 30 Mj, and 0.1 M© is greater than 2.5-10%, 14-55%, and 25-100%, respectively, 
at 100 AU distance from the star. Although such enhancements are not large, they 
may be detectable in the near future with sufficiently sensitive observations. In the fast 
wind assumption, the mass of the unseen object may be estimated by measuring the 
density jump of the spiral arm at a given position, given knowledge of the sound speed. 
Furthermore, if more than one turn of the spiral is detected, the orbital properties could 
be also estimated by fitting the density peak profile. 

In this case that more than one turn are detected, additional information can be 
inferred from the arm interval. Assuming the pattern propagation speed does not vary 
significantly along one turn, the integration of equation (3) gives the arm interval, defined 
by the distance between the outer boundary of the pattern, as following: 

^rarm = (Hv + Cs ) X — — ^ . (6) 
Vp 

From this result, one can estimate the orbital period of the companion provided that 
^rarm, ywifarm), and c^{rarm) are known from multi- wavelength observations. Alterna- 
tively, it can be rewritten as 



2;r(VH, + Cs) ^ 3/2 
(GM,)V2 



which can be used to determine the orbital distance of the companion provided that 
knowledge of the stellar mass is obtained separately. 

The application of our results to the analysis of observation data is somewhat quali- 
tative given that there are many uncertainties in the inner part of the envelopes resulting 
from the effects associated with the detailed driving of the stellar wind by stellar radi- 
ation, pulsation of the AGB star, the gas heating-cooling properties, and the accretion 
properties of the substellar object. Nevertheless our results suggest that the search for 
spiral arm patterns in the circumstellar envelope of an AGB star is desirable in order to 
make use of the full potential of the probes described in this paper to infer the presence 
and to constrain the properties of unseen objects orbiting about evolved mass-losing 
stars. 
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